Nonminimal coupling of perfect fluids to curvature 
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In this work, we consider different forms of relativistic perfect fluid Lagrangian densities, that 
yield the same gravitational field equations in General Relativity (GR). A particularly intriguing 
example is the case with couplings of the form [1 + / 2 (i?)]£ m , where R is the scalar curvature, which 
induces an extra force that depends on the form of the Lagrangian density. It has been found that, 
considering the Lagrangian density C m = p, where p is the pressure, the extra-force vanishes. We 
argue that this is not the unique choice for the matter Lagrangian density, and that more natural 
forms for C m do not imply the vanishing of the extra-force. Particular attention is paid to the 
impact on the classical equivalence between different Lagrangian descriptions of a perfect fluid. 
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I. INTRODUCTION 

Recently, in the context of f(R) modified theories of 
gravity, it was shown that a function of R— matter cou- 
pling induces a non- vanishing covariant derivative of the 
energy-momentum, V M T Ml/ ^ 0. This potentially leads 
to a deviation from geodesic motion, and consequently 
the appearance of an extra force [1]. Implications, for 
instance, for stellar equilibrium in this context have been 
studied in Ref. [2]. The equivalence with scalar-tensor 
theories with two scalar fields has been considered in Ref. 
[3] , and a viability stability criterion was also analyzed in 
Ref. [4]. This novel coupling has attracted some atten- 
tion and, actually, in a recent paper [5], this possibility 
has been applied to distinct matter contents, such as a 
massive scalar field and a dust distribution. Regarding 
the latter, it was argued that a "natural choice" for the 
matter Lagrangian density for perfect fluids is C m = p, 
based on Refs. [6, 7], where p is the pressure. This choice 
has a particularly interesting application in the analysis 
of the R— matter coupling for perfect fluids, which implies 
in the vanishing of the extra force. However, we point out 
that, despite the fact that C m = p does indeed reproduce 
the perfect fluid equation of state, it is not unique: Other 
choices include, for instance, C m = —p [7, 8], where p is 
the energy density, or C m — ~na, where n is the particle 



number density, and a is the physical free energy defined 
as a — p/n — Ts, with T being the fluid temperature 
and s the entropy per particle. Indeed, all these are on- 
shell representations of a more general Lagrangian den- 
sity, that is, obtained through back-substitution of the 
equations of motion into the related action (see Ref. [7] 
for details). Furthermore, this equivalence is established 
within the framework of GR. In this work, we address 
the issue of the Lagrangian formulation of perfect fluids 
in the context of the proposed model with a non-minimal 
coupling of the scalar curvature to matter, as depicted 
below. 

This paper is organized as follows: In Section II, we 
review the equations of motion in a curvature-matter cou- 
pling; in Section III, we show the non-uniqueness of the 
relativistic perfect matter Lagrangian densities; in Sec- 
tion IV, we analyze the perfect fluid Lagrangian descrip- 
tion with a non-minimal scalar curvature coupling; and 
in Section V, we present our conclusions. Throughout 
this work, we use the convention k = 8ttG = 1 and the 
metric signature (—1,1,1,1). 



II. EQUATION OF MOTION WITH 
CURVATURE-MATTER COUPLINGS 



The action for curvature-matter couplings, in f(R) 
modified theories of gravity [1], takes the following form 
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-f 1 (R) + [l + Xf 2 (R)]C r , 



-gd 4 x, (1) 



where fi{R) (with i — 1,2) are arbitrary functions of 
the curvature scalar R and C m is the Lagrangian density 
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corresponding to matter. 

Varying the action with respect to the metric g M „ yields 
the field equations, given by 



1 



FiR^ v - -fig^ - V„V„Fi + g^UFx = (1 

n)c 



(2) 



where we have denoted Fi(R) = fl{R), and the prime 
represents the derivative with respect to the scalar cur- 
vature. The matter energy-momentum tensor is defined 
as 



T 

± Ul> 



2 Sj^Cm) 



(3) 



Now, taking into account the generalized Bianchi iden- 
tities, one deduces the following corrected conservation 
equation 



XF 2 



A/ 2 



(4) 



If one considers the equivalence with a scalar field theory 
(with two scalar fields, <j) — R and ip — C m ) [3], it is 
clear that the non-minimal coupling between curvature 
and matter yields an exchange of energy and momentum 
between the latter and these scalar fields. 

In the following, consider the equation of state for a 
perfect fluid 



Tpv = ip + p) U^U U + pg^ , 



(5) 



where p is the energy density and p, the pressure, re- 
spectively. The four- velocity, U^, satisfies the conditions 
XJ^IJV = -1 and U^U^ = 0. 

Introducing the projection operator = g^ + U^U^, 
gives rise to non-geodesic motion governed by the follow- 
ing equation of motion for a fluid clement 



ds 



K u a u fj 



where the extra force, is given by 



P + P 



1 + A/ s 



(C m ~ P) V„R + V„p 



(6) 



h"" . (7) 



One verifies that the first term vanishes for the specific 
choice of C m = p, as noted by Ref. [5]. However, as 
emphasized in the Introduction, this is not the unique 
choice for the Lagrangian density of a perfect fluid, as we 
shall outline below. 



III. RELATIVISTIC PERFECT FLUID MATTER 
LAGRANGIAN DENSITIES 

In this section, we follow Ref. [7] closely, and review 
the Lagrangian formulation of a perfect fluid in the con- 
text of GR. The action is presented in terms of Lagrange 



multipliers along the Lagrange coordinates a A in order 
to enforce specific constraints, and is given by 

S m = J d A x [-y/=g p(n, s) + J" (tp tll + + (3 A a A )] • 

(8) 

Note that the action S m — S(g^, J M , <p, 9, s, a A , (3a) is a 
functional of the spacetime metric g^, the entropy per 
particle s, the Lagrangian coordinates a A , and spacetime 
scalars denoted by ip, 9, and (3a, where the index A takes 
the values 1, 2, 3 (see Ref. [7] for details). The physical 
interpretation of these parameters is given below. 

The vector density J M is interpreted as the flux vec- 
tor of the particle number density, and defined as J M = 
^—gnU^. The particle number density is given by 
n = \ J\/V~9' so t na * the energy density is given as a 
function p = p(\J\/^/^g, s). 

Varying the action with respect to the metric, and us- 
ing the definition given by Eq. (3), provides the stress- 
energy tensor for a perfect fluid 

= p U^U V + (n^ - p) (g^ + U"U l/ ) , (9) 



with the pressure defined as 



dp 

p = n- p. 

on 



(10) 



Note that this definition of pressure is in agreement with 
the First Law of Thermodynamics, dp = udn + nTds. 
The latter shows that the equation of state can be speci- 
fied by giving the function p(n, s), i.e., the energy density 
as a function of number density and entropy per parti- 
cle. The quantity p, = dp/dn = (p + p)/n is defined 
as the chemical potential, which is the energy per par- 
ticle required to inject a small amount of fluid into a 
fluid sample, maintaining a constant sample volume and 
a constant entropy per particle s. In addition, when im- 
posing the stress-energy tensor covariant conservation, 
i.e., T.y = 0, the perfect fluid also implies the covariant 
conservation of particle number, given by (nU^)-^ = 0. 

The variation of the action with respect to J M , ip, 9, s, 
a A and (3 a, provides the following equations of motion, 
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pU„ + tp^ + s9^ + (3 A a A = , 




-(^J"), / , = o, 

a A J» = 0. 



(11) 
(12) 

(13) 
(14) 
(15) 
(16) 



The first relationship, Eq. (11), provides the velocity- 
representation of the 4— velocity; the second equation, 
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Eq. (12), reflects the particle number conservation, i.e., 
(nU^).^ = -JL-Jt^ = 0; Eq. (13) translates the entropy 

exchange constraint; Eq. (14) provides the identification 
of T = tl JJ^ = ^-ff U after comparing it with the First 
Law of Thermodynamics; Eq. (15) reflects the constancy 
of the parameter (3 a along the fluid flow lines; and finally, 
Eq. (16) restricts the fluid 4— velocity to be directed 
along the flow lines of constant a A . 

One may now infer the physical interpretation for the 
respective parameters. The scalar field tp is interpreted 
as a potential for the chemical free energy /, and is a La- 
grange multiplier for the particle number conserva- 
tion. The scalar fields (3a are interpreted as the Lagrange 
multipliers for a A ^J^ = 0, restricting the fluid 4— velocity 
to be directed along the flow lines of constant a A . 

Note that taking into account Eq. (11), and the defi- 
nitions J p = ^J~^gnU^ and p = (p+p)/n, the action Eq. 
(8) reduces to the on-shell Lagrangian density £ m (i) = p, 
with the action given by 



(17) 



which is the form considered in Ref. [6] . One should bear 
in mind that this on-shell Lagrangian density yields the 
equations of motion (11)-(16) only if one considers that 
the pressure is functionally dependent on the previously 
considered fields ip, s, 0, (3a, ce A , and on the current 
density J^. 

Now, it was a Lagrangian density given by C m = p that 
the authors of Ref. [5] use to obtain a vanishing extra- 
force due to the non-trivial coupling of matter to the 
scalar curvature R. For concreteness, replacing C m = p 
in Eq. (7), one arrives at the general relativistic expres- 
sion 



h^V v p 
p + p 



(18) 



However, the on-shell degeneracy of the Lagrangian den- 
sities arises from adding up surface integrals to the ac- 
tion. For instance, consider the following surface inte- 
grals added to the action Eq. (8), 



-Jd 4 x(pJ^). 



Jd 4 x{0sJ^)^, 



-/ d A x{J»(3 A a A ), ll , 
so that the resulting action takes the form 



/ 



d 4 x 



-g p(n, s) - p,J(l 



e{sJ»)^-a A ((3 A Jn,L 



(19) 



Note that this action reproduces the equations of motion, 
Eqs. (11)-(16). Taking into account the latter, the action 
reduces to 



/ 



d 4 x. 



(20) 



i.e., the on-shell matter Lagrangian density takes the fol- 
lowing form C m = —p; as before, p is dependent on the 
original fields present in the action Eq. (8). This choice 
is also considered for isentropic fluids, where the entropy 
per particle is constant s = const [7, 8]. For the latter, 
the First Law of Thermodynamics indicates that isen- 
tropic fluids are described by an equation of state of the 
form a(n,T) = p(n)/n - sT [7] (see Ref. [9] for a bulk- 
brane discussion of this choice). 

For this specific choice of £ m ( 2 ) = — p the extra force 
takes the following form: 



XF 2 



1 + A/ 2 



1 



p + p 



(21) 



An interesting characteristic is that the term related to 
the specific curvature-matter coupling is independent on 
the energy-matter distribution. 

Another interesting action functional is given by the 
equation of state of the physical free energy as a function 
of the number density and the temperature, a(n,T). For 
this we follow the reasoning of Ref. [7]. For instance, 
solving Eq. (14) for s as a function of n and T (using the 
definition T = # i(U J M /| J\), and finally eliminating s from 
the action Eq. (8), yields 



J d 4 x [- 



\J\ a(n,T) + J» (cp^+p A a A )} 



(22) 



Using the definitions n = 



r=g and T = 6^J»/\J\, 



and varying the action with respect to J M , one ends up 
with the equation of motion Eq. (11). The remaining 
equations of motion are readily obtained by varying S 
with respect to ip, 9 A and (3 A . It is simple to show that 
this action also provides the perfect fluid stress-energy 
tensor. As before, one may consider the addition of the 
following surface integrals to Eq. (8) 



- J d 4 2#^),u , - / d 4 x( J»(3 A a A ) 
so that the action takes the following form 
S m = J d A x (- no) . 



jX ) 



(23) 



The matter Lagrangian density is given by £ m (3) = — na. 
The extra force in terms of this Lagrangian density yields 
the following expression: 



r = 



i 



p+p 



XF 2 



1 + A/ 2 



{na + p)V v R + V„p 



h» v . (24) 



Hence, it is clear that no immediate conclusion may 
be extracted regarding the additional force imposed by 
the non-minimal coupling of curvature to matter, given 
the different available choices for the Lagrangian density; 
moreover, one could doubt the validity of a conclusion 
that allows for different physical predictions arising from 
these apparently equivalent Lagrangian densities. 
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Thus, from the above point of view, there is no par- 
ticular reason to regard the choice of the on-shell La- 
grangian density £ m (i) — p as preferable over the others 
we have discussed above. However, this degeneracy of 
the Lagrangian density of a perfect fluid, which does not 
appear in GR, is rather intriguing and object of further 
discussion in the next section. 



IV. PERFECT FLUID LAGRANGIAN 
DESCRIPTION WITH NON-MINIMAL SCALAR 
CURVATURE COUPLING 

There is a caveat in above treatment, which can easily 
pass ignored: The discussion of the Lagrangian density- 
dependence of the extra force given by Eq. (7), and de- 
generacy thereof, implicitly admits that the equivalence 
between different on-shell Lagrangian densities holds. 
However, the latter is established in GR, and may not 
be valid in the more general model considered here. 

Clearly, one may argue that two Lagrangian densi- 
ties are equivalent if both generate the same energy- 
momentum tensor, and if variation of the corresponding 
actions yields the same equations of motion (11)-(16). As 
has been shown above, the several on-shell Lagrangian 
densities £ m (i) = p, £ m (2) = — P-, £m(3) = —na are all 
equivalent to the original, "bare" Lagrangian density, 



C m = -p(n, s) 



s6 , 



(3 A a J 



(25) 



Hence, one must attempt to retrace the derivation of the 
classical equivalence leading to these on-shell quantities. 
Clearly, if one simply includes the [1 + A/2(i?)] factor of 
Eq. (1) into action Eq. (8), that is, 

S = J d 4 x[l + \f 2 (R)} x 

[-V=g p(n, s) + J" (<p tli + 80 tli + A a%)] (26) 

then the equations of motion (11)-(16) are unaffected, as 
variation with respect to each field yields only a global 
factor [1 + A/ 2 (.R)]. 

However, the guiding principle behind the proposal 
first put forward in Ref. [1] is to allow for a non-minimal 
coupling between curvature and matter. Thus, the mod- 
ification of the perfect fluid action Eq. (8) should only 
affect the terms that show a minimal coupling between 
curvature and matter, i.e., those multiplied by y/— g. For 
this reason, the current density term, which is not cou- 
pled to curvature, should not be altered. This yields 



./a 



-V=g[l + \f 2 {R)]p{n,8) 



(27) 



The equations of motion (12), (13), (15) and (16) are 
unchanged, while Eqs. (11) and (14) read 

SS 



fj. [1 + h{R)\ Up + <p llt + s0,„ + f3 A a A 



0, 



(28) 



This results from the coupling of the variables J M and s 
with the factor [1 + f 2 (R)] (since n = \J\/y/—g). Recall- 
ing that J M = yf—gnU^, one obtains 

(29) 
(30) 



[1 + A/ 2 (#)] pUf, = <p tli + s6 tli + p A a A 



T = 



ldp 

n ds 



I + A./2OR) 



so that both the velocity representation and the temper- 
ature reflect the non-minimal coupling of curvature to 
matter. 

One may now proceed and substitute the modified 
equations of motion into action (27), in order to obtain 
the new on-shell Lagrangian density, 



S m 



(31) 



Hence, one concludes that the on-shell Lagrangian den- 
sity >C m (i) — p is also obtained in the considered scenario. 
By including extra surface integrals, a similar procedure 
(not pursued here) also yields the previously discussed 
forms C m( 2) = ~P, £m(3) = ~na. 



A. Gravitational field equations and the 
nonequivalence between on-shell and bare 
Lagrangian densities 

The above discussion confirms that one may adopt any 
particular on-shell Lagrangian density as a suitable func- 
tional for describing a perfect fluid, therefore leading to 
the issue of distinguishing between different predictions 
for the extra force. However, this is not quite correct: 
although the above Lagrangian densities C m u\ are in- 
deed obtainable from the original action, it turns out 
that they are not equivalent to the original Lagrangian 
density C m . Indeed, this equivalence demands that not 
only the equations of motion of the fields describing the 
perfect fluid remain invariant, but also that the gravita- 
tional field equations do not change. 

Recall that the terms in the field equations (2) which 
depend on C m arise from the presence of the non- minimal 
coupling [1 + A/2(i?)]. However, the formulation of a 
perfect fluid action functional includes the presence of a 
current density term, plus eventual surface integral terms 



Writing C c 
simplicity, then 



-p(n, s) , Vfj, = ip^ + s6 tl j, + , for 



Xf 2 (R)]jC c + J^ 



(32) 

one can see that only the non-minimal coupled term C c 
appears in the field equations, as variations with respect 
to g^ v of the remaining terms vanish: 

FiR^ - -hg^u - V^V.Fx + g^UFi = (1 + A/ 2 )T^ 



-2\F 2 £ C R^ + 2A(V M V, - g^U)C c F 2 , 



(33) 
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Clearly, the appropriate energy-momentum tensor is still 
obtained from C c , definitions (fO) and relations U^U^ = 
-1 and J» = ^gnU^. 

One arrives not at a paradox, but a tautology: different 
predictions for non-geodesic motion result from different 
forms of the gravitational field equations. Therefore, the 
equivalence between different on-shell Lagrangian den- 
sities £ TO (j) and the original quantity C m is broken, so 
that one can no longer freely choose between the avail- 
able forms. 

By the same token, the additional extra force is unique, 
and obtained by replacing C c = —p into Eq. (7), yielding 
expression (21), here repeated for convenience 



obtaining 



r 



XF 2 



1 + A/ : 



■V V R- 



1 



P + P 



(34) 



B. Null dust case 

Following Ref. [5] , it is interesting to analyze the gener- 
alized conservation law, given by Eq. (4) , in the case of a 
null dust matter distribution. The latter is defined as the 
particular case of a perfect fluid with vanishing pressure, 
p = 0. This is usually interpreted as expressing weakly 
interacting non-relativistic particles, with pc 2 »psi0. 
However, given the previous discussion of the functional 
description of a perfect fluid, where the pressure is not an 
independent quantity, but defined by Eq. (10), a more 
rigorous (and physically compatible) formulation corre- 
sponds to an isentropic (s = const) perfect fluid with an 
equation of state of the form p{n) = n/x, with a constant 
chemical potential p. 

The authors of Ref. [5] conclude that the extra force 
arising due to the non-minimal coupling of dust with 
the scalar curvature docs not lead to non-geodesic mo- 
tion, as it preserves parallel transport (and only changes 
the parameterization of the geodesic). However, this re- 
sult arises from the particular choice C m = p, which is 
commonly used in the framework of GR. As the previ- 
ous discussion has shown, in the context of the consid- 
ered curvature-matter coupling model, one cannot freely 
chose between available on-shell Lagrangian densities, 
since these do not lead to the same gravitational field 
equations [15]. 

Instead, inserting the component of the original La- 
grangian density that couples to the geometry, C c = —p 
into Eq. (4), and considering the energy- momentum ten- 
sor = pUf^Uv, one arrives at the following relation- 
ship 

(U»V^U V + UvVUn + UJJ^) P = 

-Tzrrbvu + UvUjwR. (35) 

1 + A/2 

Following the notation of Ref. [5] , one writes 

6 - V% + -U^p + -^L-U^R , (36) 

P 1 + A/2 



uoVuPv = -eu v 



XF 2 



1 + A/ 2 



V V R, 



(37) 



which clearly shows that parallel transport is no longer 
conserved, and one concludes that non-geodesic motion 
is also followed by pressureless dust. 



V. CONCLUSIONS 

In this work we have discussed the degeneracy of La- 
grangian densities for a perfect fluid, in the context of 
a gravity model where matter is coupled non-minimally 
with the scalar curvature. This degeneracy problem is 
well known in the context of GR, but in the discussed 
non-minimally coupled model possesses some new and 
rather surprising features, such as non-geodesical motion 
(first discussed in Ref. [1]; see Ref. [10] for a recent re- 
view); this dependency on the choice of the Lagrangian 
density was pointed out in Ref. [5]. 

We show that this degeneracy does not appear in 
the considered model, since different on-shell Lagrangian 
densities which are classically equivalent do not yield the 
same gravitational field equations. Instead, we conclude 
that only the part of the original Lagrangian density C m 
that is coupled to the geometry (via the \f^g factor) 
appears in these field equations. Hence, it follows that 
the motion of test particles is necessarily non-geodesic, if 
the non-minimal nature of the coupling between matter 
and curvature is properly accounted in the onset of fluid 
treatment. 

However, we should point out that this study only 
solves the issue of an apparent degeneracy due to the 
classical equivalence between on-shell Lagrangian densi- 
ties. This should not be seen as an exhaustive account 
of the overall problem, since it only lifts this degeneracy 
for a particular original Lagrangian density C m . One 
can take a Lagrangian density different from that of Eq. 
(8) to begin with, which also gives a full account of the 
behavior of a perfect fluid, describing both the correct 
energy-momentum tensor as well as its thermodynamics. 
As an example, if the bare Lagrangian density of Ref. 
[11] is considered (see Ref. [12] for a discussion), one ob- 
tains an extra force that is corrected by a factor linearly 
dependent on the Hclmholtz free energy. 

If one takes an initial C m that is functionally different 
from the one adopted in this study, and that still enables 
a convenient description of a perfect fluid (and is suitably 
interpreted through the use of the First Law of Thermo- 
dynamics), then one could obtain different results for the 
predicted extra force. This might lead to two different 
interpretations: One can conjecture that there must exist 
an underlying principle or symmetry that yields a unique 
Lagrangian density for a perfect fluid, so that different 
extra force predictions stem from an incomplete action 
description of it; however, one might also posit that dif- 
ferent extra forces arising from different Lagrangian den- 
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sities are physically distinguishable. If so, the model un- 
der scrutiny would serve to discriminate between differ- 
ent fluids that share the same energy-momentum tensor 
(and are thus "perfect"), but have different thermody- 
namic formulations. In the authors' opinion, this has 
not been given due attention in the literature, most likely 
because arbitrary gravitational field equations depending 
on the matter Lagrangian have not often been the object 
of scrutiny. 

In fact, to judge which matter Lagrangian density is 
the "natural" one depends, to some extent, on the inde- 
pendent variables that are considered. In this respect, 
an interesting avenue for future research would be to 
consider the non-minimal curvature-matter coupling us- 
ing velocity-potentials [6]. In this case, there are no 
constraints in the action principle, and one could com- 
pare this analysis with the one where the comoving La- 
grangian coordinates label the fluid elements and that 
exhibits constraints [13] (we refer the reader to Ref. [14] 
for a proof of the need for constraints and related issues). 



However, this is not a trivial task, as no Lagrangian is 
unique, even in the presence of the non-minimal coupling, 
since it is invariant under the addition of a divergence, 
as mentioned above. Work along these lines is presently 
underway. 
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